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ABSTRACT 

In this paper we use one of the Post-Keplerian parameters to obtain constraints on 
/(i?)-theories of gravity. Using Minkowskian Hmit, we compute the prediction of f{R)- 
theories on the first time derivative of the orbital period of a sample of binary stars, 
and we use our theoretical results to perform a comparison with the observed one. 
Selecting a sample of relativistic binary systems we estimate the parameters of an 
analytic /(i?)-gravity. We find that the theory is not ruled out if we consider only the 
double neutron star systems, and in this case we can cover the existing gap between 
the General Relativity prediction and the observed data. 
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1 INTRODUCTION 

The gravitational waves (GWs) are one of the most promis- 
ing tools to study astrophysical systems like Neutron Stars 
(NS), coalescing binary systems. Black Holes (BHs), and 
White Dwarfs (WDs). 

The observational indirect evidences of gravitational ra- 
diation were measured on the system B1913-I-16, known as 
the Hulse- Taylor binary pulsar, and confirmed in others rela- 
tivistic binary systems. The prediction of General Relativity 
(GR) on the first time derivative o f the orbital period in bi - 
nar y pulsar systems was st udied bv lHulse and Tavlod (|l975l l 
and lWeisberg et al.l (|2010l ). for which the discrepancy on ob- 
served data with respect to the prediction is ~ 1%. However, 
the observational results sho uld be explained using a dif- 
fe rent formulation of gravity ([FYeire et al]|2012l '). As shown 
in lDe Laurentis and Capozziellol ^ 201 j ). tliese systems could 
represent a good test for Extended Theories of Gravity 
(ETG). Considering a class of analytic /(i?)-theories, it is 
possible evaluate the gravitational radiated power in weak 
field limit. In this approximation we find that the energy 
radiated depends on the third derivative of the quadrupole, 
as predicted by GR, and the fourth derivative representing 
the corrective contribution to the theory. This result can be 
used to set constraints on the theory, comparing the predic- 
tion on the first time derivatives of the orbital period with 
respect to the observed one. The outline of the paper is the 
following: in sec. 2 we briefiy introduce the weak field limit 
approximation of /(_R)-theories of gravity. In sec. 3 we apply 
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the theoretical results previously obtained to binary systems 
computing the energy lost through GWs emission. In sec. 4 
we compute the first time derivative of the orbital period 
in /(i?)-theories of gravity, and we compare the theoretical 
prediction with the observed data. Finally, in sec. 5 we give 
our conclusions and remarks. 



2 F(R)-GRAVITY BACKGROUND 

The /(i?)-theories are based on corrections and enlarge- 
ments of the GR theory adding higher-order curvature in- 
variants and minimally or non-minimally coupled scalar 
fields into dynamics which come out from the effective act ion 
of quantum gravity (|Capozziello and De Laurentis! [20T1I ). 

Starting from the following field equa- 
tions in /( fl)-gravity (looking for major de- 
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i200i)S 

SnfiR) + f'{R)R - 2f{R) = |r , (2) 

the Minkowskian limit can be calculate for a class of analytic 
/(7i)-Lagrangian (i.e. Taylor expandable in term of the Ricci 
scalai0) 



(R-RoT ~fo + f^R+^R' + 



(3) 



At the first order, in term of the perturbations, the field 
equations become 



R^^ 
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where /o = ^1 , fo = -J^l and □ = ,<^''^ that is 
d'Alembert operator. Here, the Ricci tensor and scalar read 



rW = /i^.^.<^^ - D/i 



(5) 



Now, assuming that the source is localized in a finite region, 
as a consequence outside this region T^^ — and 



(9) 



To simplify the above equation the weak field limit ap- 
proximation is taken into account, i.e. the source /i^^ will be 
written as function of a single scalar variabl e t' = t — r, and 
as a consequence, it will be almost plane (|Maggiorel l2007l : 
iDe Laurentis and Capozziellci|[201 j ). 

Finally, the energy momentum tensor assumes the fol- 
lowing form 



GR 



^ V 



(10) 



To be more precise, the first term, depending on the choice 
of the constant /o, is the standard GR term, the second is 
the f{R) contribution. It is worth noticing that the order of 
derivative is increased of two degrees consistently to the fact 
that /(fl)-gravity is of fo urth-order in the metric approach 
l|De Laurentis and CapozziellQi20lll ). 
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From here it is possible calculate the energy momen- 
tum tensor of gravitation al field in /(J?)-gravity, a dopt- 
ing the definition given in iLandau and Lifshit j (|l962l ) and 
iDe Laurentis and Capozziellol (|2011^ . so that it satisfies a 
conservation law as required by the Bianchi identities: 
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Starting from ab ove equation, 

iDe Laurentis and Capozziellol (|201ll ) have shown that 
the energy momentum tensor consists of a sum of a GR 
contribution plus a term coming from /(7?)-gravity 



that in term of the perturbation h is 



(8) 



3 RADIATED ENERGY 

In order to c alculate the radiated energy of a grav itational 
waves source. IPe Laurentis and Capozziellol (|201ll ') suppose 
that hpi, can be represented by a discrete spectral represen- 
tation. The periodicity T will be proportional to the inverse 
of the difference of the pair of frequency components in the 

wave, and then, the average of must be evaluated over 

an in t erval equa l to or greater than T (|Landau and Lifshit zj 
1 19621 : lMaggior3 l2007l ). The instantaneous fiux of energy 
through a surface of area r^dQ, in the direction x is given by 



— — = r^dQx^t'^^ , 
dt 

and the average flux of energy can be written as 
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Defining the following moments of the mass-energy dis- 
tribution: 
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signature, and indicate with "," partial derivative and with " ;" 
covariant derivative with regard to gjji^; all Greek indices run from 
0, ...,3 and Latin indices run from 1, ...,3; g is the determinant. 
^ For convenience we will use / instead of f{R). All considera- 
tions are developed here in metric formalism. From now on we 
assume physical units G = c = 1. 
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and analyzing the radia t fon in terms of multipoles, 
iDe Laurentis and Capozziellol (|201ll ) found 
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iMaggiord (|2007l ) 
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Using the result in eq. (|12p and integrating over all di- 
rections they computed the total average flux of energy due 
to the tensor wave, 
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Precisely, for /o' and /o — >■ I, eq. p7[l becomes 



/ Gnic 
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(20) 



where a is the semi-major-axis, and e is the orbital eccentric- 
ity. We obtain the following relations for the time derivatives 
of the quadrupole 



Qii = Hi sin 2?/i(ecos ?/) -f 1) (3e cos ?/) -I- 4), 

Q22 = -Hi(8cos7/) -I- e(3cos27/) -f 5)) X 
X sin?/'(ecos7/; + 1)^, 

Q12 = -Hi(ecos7/) -f 1)^ X 

X (5e cos xf) + 3e cos Stl) + 8 cos 2tp) , 



(21) 
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which is t he well-known result of GR l|Landau and Lifshit j 
1 19621 : IWeinbcra ,197 2). An important remark is related to 
the absence of monopole and dipole terms in our consid- 
erations. In our case, all the calculations are performed in 
the Jordan frame so/(JfI)-gravity results as a mere exten- 
sion of GR being f(R) = R, so any dipole terms is null 
(as shown in IWilll (119931 ) in table 10.2). In order to put in 
evidence such contributions, we have to pass in the Einstein 
frame where the additional degrees of freedom of gravita- 
tional field can be recast in term of scalar fields. I n this case , 
monopole and dipole t erms explicitly come out (Will 1993; 
iNaef and Jetzed I2OIII : |d amour and Esposito-Farcsc 1 9961 ). 
The two approaches are conformally equivalent but in the 
Einstein frame monopole and dipole terms can come out 
(see, e.g. ICapozziello and De Laurentid (|201ll ')'). 



4 APPLICATION TO PULSAR BINARY 
SYSTEMS 

Now, our goal is to use a sample of binary pulsar systems to 
fix bounds on f{R) parameters. To do this, we assume that 
the motion is Keplerian and the orbit is in the (x, y)-plane. 
We define nip as the pulsar mass, rric as the companion mass. 



and fj. — 



TTlc + nip 
quadrupole matrix is 



as the reduced mass. In {x, ?/)-plane the 



Qij = fJ-r 



cos tp sin ip cos 
sin tp cos ip sin^ tp 
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where i,j are the indexes in the orbital plane, r is the 
equation of the elliptic Keplerian orbit and ip is eccentric 
anomaly, and both of them are time dependent. 

To compute the radiated power we need the third 
and fourth derivatives of quadrupole, so we must compute 
the time derivatives using the following relation given in 



Q II = [l5e cos4'i/' -f 50ecos3V'+ 

+ {\2e + 32) cos 2^) + 6e cos V - 3e^] , (24) 

Q22 = — [l5e^ cos4V' 4- 50ecos3'(/'+ 

-f (24e^ -f 32) cos 2^ -f 14e cos i/> - 7e^] , (25) 



where 
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+ (33e^ + 32) cos V' + 30e] , 
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Now, from eq. (|17|) . we can perform the time average of 
the radiated power writing 



(it T 



dt 



T 



d^j dE{ip) 
ijj dt 



(27) 



and finally we get the first time derivative of the orbital 
period 

T = 3 f ryi fiGi{m, + mp)^ 
20\2nJ c5(l-e2)i 



fo {S7e* + 292e^ + 96) 



2(1 + e2)'' 



X (891e** + 280166" + 82736e'' + 43520e^ + 3072)] . 

(28) 

In the next section we will go on to constraint the f{R)- 
theories estimating /q' from comparison between the theo- 
retical predictions of Tt and the observed one. 
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4.1 Comparing theory prediction with data 

It is well known that in the relativistic binary pul- 
sar systems there is a loss of energy due to GWs 
emission. This energy loss, provided by GR, has been 
confirmed by the timing data analysis on the well 
known binary pulsar B1913-I-16 (jHulse and Tavloil 1 19751 : 
IWeisberg et alj |201C( ) . We also must note that the sys- 
tems like B1913-I- 16 are optimal tools to constrain theo- 
ries of gravitation (|Damour and Esposito-Fareselll99^ using 
the Post-Keplerian parameters. For sake of convenience we 
choose the observed orbital period derivative Tb, because it is 
one of the best observed Post-Keplerian parameters. More- 
over we know, according to GR theory, that it is related to 
the foreseen orbital decay due to quadrupole gravitational 
radiation emitted by the binary systems. 

As shown in sec. [l] it is possible rewrite the first deriva- 
tive of the orbital period in /(_R)-theories of gravity. In prin- 
ciple, if we know exactly which Lagrangian we have to use to 
describe those type of systems, then we can predict the en- 
ergy loss through GWs radiation. Here, we want to make the 
inverse process, to get an estimation of the second deriva- 
tive /o' imposing the strong hypothesis that the difference 
between the observed binary period variation {Tb^^^ ± 5) 
and the one obtained by the relativistic theory of gravita- 
tion, ATb = TbQ^^ — Tgr, is fully justified imposing that: 

^fcobs - Tgr - /o'ri,^(_R) = 0, (29) 

'^bois ±3 ~ Tgr - /(','^jT6^(„, = 0, (30) 

and propagating the experimental error, ±5, on the first 
derivative of the observed orbital period Tb^^^ , into an un- 
certainty on second derivative of gravitational theory, fg^^ . 
What we want to emphasize is that, where GR is not able 
to fully explain the loss of energy by emission of GWs radia- 
tion then, the additive contribution of an ETGs can provide 
a way to fill the gap between theory and observations. We 
also have substracted the external contributions to the or- 
bital decay as galactic or Shklovskii acceleration when those 
values are available in literature. Solving the eqs. (|29|l and 
(|30p for /o' and /q'^^ we get an estimation of /q' and its 
upper and lower limits corresponding respectively to =1=5. In 
this way AT;, is fully explained through the orbital period 
correction due to the ETGs Tt^^^j . So we get: 

fo=7p^, (31) 

/oL. - (32) 

where ATb^, = fbo,, ±S- Tgr. 

Thus, among the various binary stars catalogues avail- 
able in literature, we choose a sample of Observed Rela- 
tivistic Binary Pulsars (ORBP) such that the binary pe- 
riod r^Q^^ , the observed orbital period variation Tb^^^ , the 
computed orbital period variation from general relativistic 
theory Tgr, the orbital eccentricity e, the masses of the com- 
ponents rUp and rric, are known with a fairly good precision. 
For each system we have chosen, all previous parameters and 
their references are reported in Tab. [1] where we show: the 
J-Name of the binary pulsar system, the observed orbital bi- 



nary period TbQ^^ in days, the orbital projected semi-major 
axis asin{i) in light second, the orbital eccentricity e, the 
observed time variation of the orbital period Tb^^^ , the pre- 
dicted one Tgr (according to the GR theory), the experi- 
mental error ±5 on Tb^^^ and the masses rUp and ttIc of the 
binary system components in solar mass unit. 

Furthermore, in Tab. [2]we reported: the J-Name of the 
systems, the difference ATgr between Tb^^^ and Tgr (equal 
to the correction — /o'Tb^j^^ ), the correction Tbj^^^^ , the cor- 
responding /o' solution of (|29|) shown in pif) . the interval 
centered on /q' and finally, the interval centered on /q' and 

computed from the difference: — ^^^-^ — —, where /o'_^^, are 
the solutions of (|29[) shown in (|3ip taking in to account the 
experimental errors ±5 on the observed orbital period vari- 
ation Tto^^ . 

Now, in Fig[l] we report representative results of our nu- 
merical analysis on the sample of binary pulsars we choose. 
In both panels we use the following notation: the black line 
shows the behavior of the first derivative of the orbital bi- 
nary period for the /(_R)-theories of gravity as computed in 
eq. H28p : the blue line represents the observed orbital period 
variation Tb^^^ ; the red lines give the error band determined 
by the experimental errors ±5; and finally the green line is 
representative of the Tgr orbital period variation computed 
from the GR. For the binary pulsar system J2129 -I- 1210C, 
Fig[T^, the orbital period variation Tgr, computed from the 
GR, is included in experimental error band ±5, so as the 
observed orbital period variation Tt^^, . Moreover, it is pos- 
sible to see from Fig[T^ the GR value of Tgr is recovered for 
/o' = (green square), whilst to justify the difference ATgr 
between Tb^^^ and Tgr we have, from the solution of (|3ip . 
the values shown in Fig[T] for /o'_j_j (red square) and for /g' 
(blue square). In the panel (b) of the Fig[T] there is reported 
for J0751 -|- 1807 the same situation. In this case the Tgr is 
out of the error band determined by the experimental errors 
±(5. It is again possible to see for /y' = that the GR value 
of Tgr is recovered, but in this case the /q' values are order 
of magnitude greater than the one of the well behaved case 
of J2129 + 1210C. 

In Fig.[2]there are shown, for sake of convenience, in log- 
arithmic scale, the absolute values of /o' reported in Tab[2] 

versus the ratio X'^'" ■ We must note that for the first six bi- 

naries in tables, the ETGs are not ruled out 0.04 ^ /q' ^ 38. 

For those systems we get 0.5 ^ ^ 1.5, the difference 

between Tgr and Tb^^^ can be explained adding a new con- 
tribution from the theory of gravity. Instead for most of bi- 
naries we have /o' values that can surely rule out the theory, 
since taking account of the weak field assumption we obtain 
38 /o' ^ 4 X 10^. From this last values to the first ones, 
there is a jump of about four up to five order of magnitude 
on /('/. The origin of these strong discrepancies, perhaps, is 
due to the extreme assumption we made, to justify the dif- 
ference between the observed Tbo^^ and the predicted Tgr 
using the ETGs. 



5 DISCUSSION AND REMARKS 

We want point out in this preliminary work that, where 
the GR theory is not enough to explain the gap between 



Testing f{R)-theories using the first time derivative of the orbital period of the binary pulsars 5 



the data and the theoretical estimation of the orbital de- 
cay, there is the possibility to extend the GR theory with 
a generic f{R)- theory to cover the gap. Here, we sim- 
ply verify that this possibility exists, but there is need to 
compute the Post-Keplerian parameters in the /(7i)-theory 
to estimate correctly the masses of the binary systems to 
constraint correctly the analytic parameters the ETGs. In 
post-Minkowskian limit of analytic /(_R)-gravity models, the 
quadrupole-radiation depends on the masses of the two bod- 
ies, on the orbital parameters and on the analytic parame- 
ters of the /(i?)-theory as the coefflcients /q and /o' of the 
Taylor expansion. A first result we present is the analytical 
solution of the quadrupole radiation rate in which it is pos- 
sible separate the GR contribution and the one due to the 
/(i?) -gravity. We should note that the correction depends on 
the eccentricity of the orbit and on the orbital period of the 
binary system, and specifically, the radiation rate is a func- 
tion of /o and /o . According to eq. (P5)) ,we have selected 
a sample of relativistic binary systems for which the first 
derivative of the orbital period is observed, we have com- 
puted the theoretical quadrupole radiation rate, and finally 
we have compared it to binary system observations. From 
Tab.m it is seen that the first five systems have masses de- 
termined in a manner quite reliable, while for the remaining 
sample, masses are estimated by requiring that the mass of 
the pulsar is 1.4Mq and, assuming for the orbital inclination 
one of the usual statistical values (i = 60° or i = 90° ) , and 
from here comes then the estimate of the mass of the com- 
panion star. So a primary cause of major discrepancies, not 
only for the ETGs, but also for the GR theory, between the 
variation of the observed orbital period and the predicted 
effect of emission of gravitational waves, could be a mistake 
in the estimation of the masses of the system. In addition, 
other causes may be attributable to the evolutionary state of 
the system, which, for instance, if it does not consist of two 
neutron stars may transfer mass from companion to the neu- 
tron star. In our sample, there are only five double NS that 
can be used to test GR and ETGs. Taking into account of 
the strong hypothesis we made, the extended theory correc- 
tion to Tg r can also include the gala c tic acceleration term 
correc tion (jPamour and TavloJ(|l99lf l.lD amour and Tavlod 
Here, we give a preliminary result about the energy 
loss from binary systems and we show that, when the nature 
of the binary systems can exclude energy losses due to trade 
or loss of matter, then, we can explain the gap between the 
first time derivative of the observed orbital period and the 
theoretical one predicted by GR, using an analytical f{R)- 
theory of gravity. In conclusion, to improve the estimation 
of the /(_R)-coefficients, we need: to consider the hydrody- 
namic effects due to the transfer of the matter in the binary 
system, in order to analyze different systems from double 
NS; and to improve the estimations of the mass of the stars 
in the binary systems without prior on pulsar mass and or- 
bital inclination. 
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Table 1. Data for binary Relativistic pulsars: in the order we reported the J-Name of the binary pulsar system , the orbital binary period Tj, in days, the orbital projected semi-major 
axis a{sini) in light second, the orbital eccentricity e, the Observed orbital period variation Tobs, the predicted Tqh according to the General Relativity theory, the experimental error 
±5 on Toba a-iid the masses nip and rUc of the binary system components in solar mass unit. 



Name 


Tb 
fdavs) 


6t 

(Isec) 


i 

(degrees) 




Th 

Oba, 


(10-12) 


±5 

(10-12) 


TTlp 


TYlc 


J2129+1210C 


0.335282049 


2.51845 




0.681395 


-3.96 


-3.94 


0.05 


1.358 


1.354 


J1915+1606 


0.322997449 


2.341782 


54.12° 


0.6171334 


-2.423 


-2.403 


0.001 


1.4398 


1.3886 


J0737-3039A 


0.102251562 


1.415032 


88.69° 


0.0877775 


-1.252 


-1.248 


0.017 


1.3381 


1.2489 


J1141-6545 


0.197650959 


1.858922 


73° 


0.171884 


-0.403 


-0.387 


0.025 


1.27 


1.02 


J1537-I-1155 


0.420737299 


3.7294626 


78.4° 


0.2736767 


-0.138 


-0.192 


0.0001 


1.3332 


1.3452 


J1738-I-0333 


0.3547907399 


0.343429 


32.6° 


3.4e-7 


-0.017 


-0.0277 


0.0031 


1.46 


0.181 


J0751-I-1807 


0.263144267 


0.3966127 


65.8° 


0.00000071 


-0.031 


-0.017 


0.009 


1.7 


0.67 


J0024-7204J 


0.120664938 


0.0404021 


60° 





-0.55 


-0.03 


0.13 


1.4 


0.024 


J1701-3006B 


0.144545417 


0.2527565 


84.7° 





-5.12 


-0.09 


0.062 


1.4 


0.14 


J2051-0827 


0.099110251 


0.045052 


30° 





-15.5 


-0.03 


0.8 


1.4 


0.027 


J1909-3744 


1.533449475 


1.8979910 


86.4° 


1.302E-07 


-0.55 


-0.003 


0.03 


1.57 


0.212 


J1518+4904 


8.634005096 


20.044002 


< 47° 


0.24948451 


0.24 


-0.001 


0.22 


1.56 


1.05 


J1959+2048 


0.381966607 


0.0892253 


65° 





14.7 


-0.003 


0.8 


1.4 


0.022 


J2145-0750 


6.83893 


10.164108 




0.0000193 


0.4 


-0.0005 


0.3 


1.4 


0.5 


J0437-4715 


5.74104646 


3.36669708 


137.58° 


0.00001918 


0.159 


-0.0004 


0.283 


1.76 


0.254 


J0045-7319 


51.169451 


174.2576 


44° 


0.807949 


-3.03E-I-5 


-0.02242 


9E-I-3 


1.4 


8.8 


J2019+2425 


76.51163479 


38.7676297 


63° 


0.00011109 


-30.0 


-0.000006 


60.0 


1.33 


0.35 


J1623-2631 


191.44281 


64.80946 


40° 


0.02531545 


400.0 


-0.000003 


600.0 


1.3 


0.8 
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Table 2. Upper Limits of /q' correction to Tqh of binary relativistic pulsars assuming that all the loss of energy is caused by Gravitational 
Wave emission. Wc reported the J-Name of the system, the difference ATgh between T^q^^ and Tgr equal to the correction — /q'To^j^j , 

the correction ^b^j^jj i the corresponding /q' solution of II29I I shown in ( I31I I. the interval centered on /q' and computed from the difference 

— — — , where /q'^^ , are the solutions of ( I29|l shown in ( I31|l taking account of the experimental errors ±5 on the observed orbital 
period variation Tt^j^^ . 



Name 






/(')' 




J2129+1210C 


-2.17E-14 


6 


.OlE- 


■13 


3.61E-02 


8.32E-02 


J1915+1606 


-2.04E-14 


2, 


.lOE- 


13 


9.74E-02 


4.77E-03 


J0737-3039A 


-4.23E-15 


1, 


.86E- 


■14 


2.28E-01 


9.15E-02 


J1141-6545 


-1.65E-14 


3 


.88E- 


■15 


4.25E+00 


6.44E+00 


J1537+1155 


5.39E-14 


1 


.42E- 


■15 


-3.79E+01 


7.03E-02 


J1738+0333 


-1.56E-15 


1 


.06E- 


■16 


-1.47E+01 


2.92E+01 


J0751+1807 


1.41E-13 


8, 


.98E- 


■16 


-15.7E+01 


1.002E+01 


J0024-7204J 


-5.22E-13 


3, 


.13E- 


■16 


1.67E+03 


4.15E+02 


J1701-3006B 


-5.03E-12 


8 


.81E- 


■16 


5.71E+03 


7.04E^f01 


J2051-0827 


-1.55E-11 


4 


.77E- 


■16 


3.24E+04 


1.68E+03 


J1909-3744 


-5.47E-13 


2, 


.62E- 


■18 


2.09E+05 


1.14E+04 


J1518+4904 


2.41E-13 


3, 


.42E- 


■19 


-7.05E+05 


6.43E+03 


J1959+2048 


1.47E-11 


1 


.07E- 


■17 


-1.38E+06 


7.51E^f04 


J2145-0750 


4.01E-13 


1 


.OOE- 


■19 


-4.00E+06 


2.99E+06 


J0437-4715 


1.59E-13 


1, 


.04E- 


■19 


-1.57E+06 


2.73E+06 


J0045-7319 


3.02E-07 


1, 


.llE- 


■16 


2.74E+9 


8.13E+07 


J2019+2425 


-3.00E-11 


1 


.llE- 


■22 


2.71E+11 


5.41E-f 11 


J1623-2631 


4.00E-10 


2 


.02E- 


■23 


-1.98E-f 13 


2.97E+13 
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J2129+1210C 




(a) 



J0751+1807 




(b) 

Figure 1. We report representative results of our numerical analyses on ttie sample of binary pulsars we have selected. In both figures 
we use the following notation: the black line shows the behavior of the first derivative of the orbital binary period for the /(ij)-thcory of 
gravitation as computed in cq. I I28I I: the blue line represents the observed orbital period variation Ti,^^^; the rod lines give the error band 
determined by the experimental errors ite; and finally the green line is representative of the Tqh orbital period variation computed from 
the GR. In the panel (a) for the system J2129+1210 C the Tgr is included in the error band determined by the experimental errors 
±5, so as . We point out the GR value of Tgr is recovered for f"{r) = (green square), while to justify the difference between 

''obs '^G-R we show the value of /q' (blue square) and its error band /g'^^ (red square) as computed in cqs. II31I I and I I32II . In the 
last panel (b) there arc reported for J0751+1807 the same data but in this case the Tgr is OUT of the error band determined by the 
experimental errors ±5. It is possible to see for /q' = that the GR value of Tqr is recovered, but in this case the /q' values arc much 
greater than the previous ones. 
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Figure 2. In figure there arc shown, for sake of convenience, in logaritmic scale, the absolute values of /q' reported in Tab. [2] versus 



the ratio ^ . O^^ , We must note that for five binaries the ETGs v^e are probing is not ruled out 0.04 5^ fl^ ^Si 38, for those systems the 
Tgr 

difference between Tgr and Tbobs indeed we get 0.5 ^ ^ 1.5. Instead for most of binaries we have /q' values that can 

surely rule out the theory, since taking account of the weak field assumption we obtain 38 /q ^ 4 X lO'^. From this last values to the 
first ones, there is a jump of about four up to five order of magnitude on /g. 



